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ABSTRACT 

We study the accretion of hot dark matter onto moving cosmic string loops, using 
an adaptation of the Zeldovich approximation to HDM. We show that a large nnmber of 
nonlinear objects of mass greater than 10^^ Mq, which could be the hosts of high redshift 
quasars, are formed by a redshift of 2 : = 4. 
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1 Introduction 


Topological defect theories provide an alternative to inflation for explaining the origin of the 
primordial fluctuations which have grown into present-day structures through gravitational 
instability. Cosmic strings are one-dimensional topological defects possibly formed in a phase 
transition in the early Universe. 

The observed presence of massive nonlinear structures at high redshifts of 3 to 4 has provided 
stringent constraints on models of structure formation containing hot dark matter (HDM), a 
candidate for which is a massive neutrino. This is due to the large thermal velocities of the 
HDM particles, which prevent the growth of perturbations on scales smaller than the free- 
streaming length. For inflationary models with adiabatic density flucUrations, recent data 
on the abundance of damped Lyman alpha absorption systems (DLAS)EJ and on the quasar 
abundancell has restricted the fraction of HDM to less than about 30% of the total dark matter 
presentB. The scenario of structure formation with cosmic strings is, however, viable even if 
the dark matter in the universe is hot, because cosmic strinra, which provide the seeds for 
the density perturbations, survive the neutrino free-streamingo B . In this article we present 
work onthe constraints imposed on the cosmic string theory by the abundance of high redshift 
quasars O. 

Searches for high redshift quasars have been going on for some timeo. The quasar luminos¬ 
ity function is observed to rise sharply as a function of redshift z until z ~ 2.5. According to 
recent results from the Palomar grism survey by Schmidt et al. (1995)0 , it peaks in the redshift 
interval ze [1.7, 2.7] and declines at higher redshifts. Quasars (QSO) are extremely luminous, 
and it is generally assumed that they are powered by accretion onto black holes. It is possible 
to estimate the mass of the host galaxy of the quasar as a function of its luminosity, assuming 
that the quasar luminosity corresponds to the Eddington luminosity of the black hole. For a 
quasar of absolute blue magnitude Mb = —26 and lifetime of tg = lO^yrs, the host galaxy 
mass can be estimated as0 

Mg = ci1012Mq, 


(1) 
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where ci is a constant which contains the uncertainties in relating blue magnitude to bolometric 
magnitude of quasars, in the baryon fraction of the Universe and in the fraction of baryons in 
the host galaxy able to form the compact central object (taken to be 10“^). The best estimate 
for Cl is about 1. Models of structure formation have to pass the test of producing enough early 
objects of sufficiently large mass to host the observed quasars. 

Besides the quasars themselves, absorption lines in their spectra due to intervening matter 
can be used to quantify the amount of matter in nonlinear structures at high redshifts. Based 
on the number density of absorption lines per frequency interval and on the column density 
calculated from individual absorption lines, the fraction of Q in bound neutral gas (denoted by 
Ug) can be estimated. For high-column density systems, the damped Ly-a systems (DBAS), 
recent observational results are that 

Ug(z) > 10-3 (2) 

for ze [1,3]. In the most recent results by Storrie-Lombardi et al. (1996)0 there is evidence 
for a flattening of Ug at z ~ 2 and a possible turnover at z ^ 3 . The corresponding value for 
ft in bound matter is larger by a factor of , where fb is the fraction of bound matter which 
is baryonic, which is about 10 % in a flat universe. 

In the next section, we give a brief review of the cosmic string scenario of structure formation. 
Then we study the accretion of hot dark matter onto cosmic string loops, which seed large 
amplitude local density contrasts already at early times. We use the results to compute the 
number density of high redshift objects as a function of a parameter v which determines the 
number density of loops in the scaling solution (see Section 2). We demonstrate that for 
realistic values of i/, the number of massive nonlinear objects at redshifts < 4 satisfies the 
recent observational constraint of quasar abundances (see Figure 1). We also comment on the 
implications of (2). We consider a spatially flat Universe containing HDM and baryons. Units 
where c = 1, a Hubble constant of H = 50 /isokms-^Mpc-^ and a redshift at equal matter and 
radiation of Zeq = 5750 are used. 

2 Cosmic Strings and Structure Formation 

Cosmic strings® are one-dimensional topological defects which might have been formed in 
a phase transition in the very early Universe. They are produced for the same reasons as 
ordinary defects in condensed matter systems, such as vortex lines in superfluid helium, but 
their subsequent dynamics is governed by relativistic equations. In the simplest model, the 
symmetry breaking occuring in the phase transition is achieved via a two-component scalar 
Higgs field (j) = {(j>i,(j) 2 ) with symmetry-breaking potential where f] is the 

energy scale of the phase transition, and A is the coupling constant. Topological considerations 
require cosmic strings to have no ends, they are either formed as infinitely long strings or as 
closed loops. These strings are very thin lines of trapped energy, the energy being the false 
vacuum energy V{(j) = 0) which is trapped inside the defects after the rest of the Universe has 
evolved to the true ground state | < ^ > | = 77 with V{(j)) = 0. The mass per unit length of the 
strings is given by the symmetry breaking scale rj , ^ ^ rf , and it is the only free parameter 
of the model (in practice uncertainties due to the complicated evolution of cosmic strings after 
the phase transition introduce extra phenomenological parameters, which are not fundamental, 
however). The energy in the defects acts gravitationally on surrounding matter and radiation, 
and thus provides the origin for the density perturbations, which evolve into today’s large-scale 
structure. 

Most important for structure formation are the infinite strings. They are straight over 
distances increasing with time of one horizon distance ~ t, and form an approximate 
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random walk on distances larger than that. In this way they can cause structures to be formed 
on the large scales observed today. The model makes predictions consistent with large-scale 
structure observationso 113 and CMB anisotropy measurements til. An intriguing fact is that 
the normalization of the free parameter ^ of the model to these observations, which gives 

« 10 “®, requires the phase transition giving rise to the defects to take place at a scale of 
77 ~ lO^^GeV . This is just the scale of grand unification where the three coupling constants 
of the strong, weak and electromagnetic interactions meet, suggesting the appearance of new 
physics. 

The network of cosmic strings is formed in a phase transition in the early universe about 
10“^® sec after the Big Bang, long before the times relevant for structure formation. Therefore 
the string network has to be evolved over many orders of magnitude in time until teq (and subse¬ 
quently) . This is very complicated numerically, and would be impossible to handle analyticallv, 
if it were not for the fact that the network of strings quickly evolves towards a scaling solutionO 
where the energy density in long strings remains a constant fraction of the total background 
energy density. This is achieved by intercommutations and self-intersections of the strings lead¬ 
ing to the production of small loops, which then decay by emitting gravitational radiation. In 
this way, some of the energy input into the string network coming from the stretching of the 
strings due to the expansion of the universe is transferred to the background. According to the 
scaling solution, the string network looks the same at all times when all distances are scaled by 
the Hubble radius and this allows the extrapolation over such long time intervals. The 

scaling solution can be pictured as having a fixed number M of long strings per Hubble volume 
at any given time, and a distribution of small loops. 

Some of the parameters characterizing the string network besides the fundamental one /i are 
the long string velocities Vs and initial loop velocities Vi, the average number M of strings per 
Hubble volume, the amount of small-scale structure on the long strings quantified by (/j, — T), 
where T is the tension of the string, and the constants a, (3 and v appearing in the next three 
equations. Small-scale structure is produced mainly when string segments self-interect and 
split off loops, a process which is characterized by the loop production rate and the size of the 
produced loops. According to the scaling solution, loops are formed with a radius which is a 
constant fraction of the horizon, 

Rf{t) = at , (3) 

and length 

I = PR. (4) 

The number of loops per unit physical volume present at time t with lengths in the interval 
between I and I + dl is 

n{l,t)dt = . (5) 


Since loops decay by emitting gravitational radiation, there is a lower cutoff value of I for the 
distribution (5) given by Imin G^t, G being Newton’s constant. Below Imin, the distribution 
niLt) becomes constant. 

Numerical simulations® indicate that a < lO’^, /3 - 10, and M ~ 10. From these values 
it follows that - unless v is extremely large - most of the mass of the string network resides in 
long strings (where long strings are defined operationally as strings which are not loops with 
radius smaller than the Hubble radius). 

T ong strings accrete matter in the form of wakes behind them as they move through 
space® 10 . Spacetime around a long straight cosmic string can be pictured as locally 
flat, but with a deficit anglefij of SnGfx . Therefore a string moving relativistically with ve¬ 
locity Vs imparts velocity perturbations to surrounding matter towards the plane swept out by 
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the string, of magnitude 

u = inGfijsVs (6) 

creating overdensities in the form of a wake behind the string. If small-scale structure is present 
on the string, there is in addition a Newtonian force towards the string, proportional to G{fi—T), 
the strings move more slowly and accrete matter rather in the form of filaments than wakes. 

For HDM, the first nonlinearities about wakes resulting from strings without any small- 
scale structure form only at late times, at a redshift oflc^ about 1 for Gfi ~ 10“®, which is 
obtained when normalizing the model to COBE. Before this redshift, no nonlinearities form as 
a consequence of accretion onto a single uniform wake. Thus, in the cosmic string and hot dark 
matter theory, a different mechanism is required in order to explain the origin of high redshift 
objects. Possible mechanisms related to wakes are early structure formation at the cro s sing 
sites of different wakesE-3, small-sc^ structure of the strings giving rise to wakesll^ 113 , 
and inhomogeneities inside of wakeslls. In this article, however, we will explore a different 
mechanism, namely the accretion of hot dark matter onto loops. 

In earlier workO, the accretion of hot dark matter onto static cosmic string loops was 
studied. It was found that in spite of free streaming, the nonlinear structure seeded by a point 
mass grows from inside out, and that the first nonlinearities form early on (accretion onto 
string filaments proceeds similarly]13). Since loop accretion leads to nonlinear structures at 
high redshift, we will now investigate this mechanism in detail to see whether enough high 
redshift massive objects to statisfy the QSO constraints and (2) form. 

3 Nonlinear mass accreted by a single static loop 

We will use the modified Zeldovich approximationlll’ 0 to study the accretion of hot dark 
matter onto moving string loops. The Zeldovich approximation is a first order Lagrangian 
perturbation theory technique in which the time evolution of the comoving displacement ip 
of a dark matter particle from the location of the seed perturbation is studied. The physical 
distance of a dark matter particle from the center of the cosmic string loop is written as 

h{q,t) = a{t)[q - ip{q,t)] . 

For a point-like seed mass of magnitude m (the string loop in our case) located at the comoving 
position q' = 0, the equation for ip is 

This equation describes how as a consequence of the seed mass, the motion of the dark 
matter particles away from the seed (driven by the expansion of the Universe) is gradually 
slowed down. If the seed perturbation is created at time ti and the dark matter is cold, then 
the appropriate initial conditions for ip are 

iP{q,ti)=ip{q,ti) = 0. (8) 

As formulated above, the Zel’dovich approximation only works for cold dark matter, particles 
with negligible thermal velocities. For HDM, free-streaming prevents the growth of perturba¬ 
tions on scales q below the free-streaming length Xj(t), which is the mean comoving distance 
traveled by neutrinos in one expansion time 


= v{t)z{t)t, 
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where v{t) ~ z{t) is the hot dark matter velocity. The free- streaming length decreases with 
time (in comoving coordinates) as so that a scale q which is initially affected by free- 

streaming, q < Xj(ti) becomes equal to Xj at some later time ts{q) , and then remains above 
it and unaffected by free-streaming at later times. 

A correct treatment of the effect of free-streaming requires the solution of the collisionless 
Boltzmann equation for the neutrino phase-space densityB. But a simple modification of the 
Zeldovich approximation for CDM gives the same results as the full treatment . For scales 
q < Xj{ti), Eq. 0 is only integrated from the time ts{q) > U onwards instead of from ti , 
i.e. growth only starts at ts(<?) , and the initial conditions are modified to tp(q,ts{q)) = 0 and 
'>p{q,ts{q)) = 0 . For scales where q > Xj{ti) initially the CDM formalism applies. 

We can now define the mass that has gone nonlinear about a seed perturbation as the rest 
mass inside of the shell which is “turning around”, i.e. for which 

h{q,t) = 0. 

In the case of CDM, this yields a nonlinear mass of 

2 f 

McDAiit) =-m\—\ , 

and for loops where the accretion is affected by free-streaming, 

8 ( t 

with Meq = 2vlgteq/{9G) , Veq is the HDM-particle velocity at teq- From Eq. || and we 
can determine, for any time t, a mass M'{t) such that for all masses M > M'{t) accreted by 
any loop so far, accretion has not been affected by free-streaming, and Eq. ^ is applicable. For 
M < expression [l^ is valid. 


(9) 


( 10 ) 


4 Number density of quasar host galaxies 


We can now compute the number density nc{> Mo^t) in nonlinear objects heavier than Mq 
at high redshift z{t) in the cosmic string and hot dark matter model (in an analogous way the 
fraction of the critical density in such objects can be calculated). By considering only the 
accretion onto strin g lo ops we will be underestimating these quantities. 

It can be shownciJ that in an HDM model string loops accrete matter independently, at 
least before the large-scale structure in wakes turns nonlinear at a redshift of about 1, i.e., later 
than the times of interest in this paper. Hence, the number density n(Z,t) of loops of length I 
given by (H) can be combined with the mass M{1) accreted by an individual loop to give the 
mass function n{M,t). Here, n{M,t)dM is the number density of objects with mass in the 
interval between M and M -|- dM at time t. This in turn determines the comoving density 
in objects of mass M > Mq- Since the functional form of M{m) and hence M{1) changes at 
M = M', the functional forms of the mass function n{M) will be different above and below 
M'. Thus 



' M'(t) 

p 

Map) 

p 

nG{> MG,t) = z ^{t) 

/ dMnniM) -\- 

/ dMnc’iM) 


_Mc 

M'(t) 
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Similarly, the fraction of the critical density in objects of mass greater than Mi, is 


^nl{t) = 


M'(t) M2(t) 

/ dMnH{M)M+ [ dMnc{M)M 


Ml 


M'(t) 


67rGf^ 


( 12 ) 


Note the upper mass cutoff M 2 {t) introduced. It is imposed by our approximation of loops as 
point masses. For masses greater than M 2 (t), this approximation breaks down, the responsible 
loops would have too large a radius to be able to accrete matter effectively. M 2 {t) can be 
estimated by demanding that the mass M (t) accreted onto a loop exceed the mass in a sphere 
of radius equal to the loop radius at ti . The integral in Eq. H is dominated at M'{t) , and can 
be approximated as 


nG{> Met) 


1 






-3 


it)- 


(13) 


In Figure 1, the comoving number density of objects massive enough to be able to host 
quasars, nG{> is plotted for G^jl = 10 “®, Veq = 0 . 1 , and the values of the parameters 

a = 10“^ and (3 = 10. The upper curve is for static loops, the lower one for initial loop velocities 
of Vi = 0.25. The effect of loop motion is to decrease the density of these massive objects, as 
discussed in the next paragraph. We have to compare our results for nGi> Mc^t) with the 
number density of host galaxies, , inferred from the observed quasar abundance nq. Since 
we assumed a finite quasar lifetime tq, the number of host galaxies is larger than nq by a 
factor of tii/tq, where tn is the Hubble time at time t, 


obs tn 
Uq = —nq . 

tQ 


(14) 


This is also plotted in Fig. 1 . ra r—i 

The formalism presented above can be adapted to moving instead of static loopsQ’ c3 . 
Loop motion changes the geometry of the accreted object, it becomes more elongated and 
the transverse scale going nonlinear becomes smaller than for loops at rest. Because of this, 
free-streaming of the neutrinos can hinder the growth of perturbations for a longer time. The 
smallest mass accreted at time t which has been unaffected by free-streaming, M'(^is increased 
compared to the case of static loops, thereby decreasing nait) according to Eq. 

The result for the fraction of nonlinear mass accreted by moving loops is 


^niit) ~ 10 '^na-2iGfi)lv^ ^z{t) 
where {Gfj,)^ is defined as (G^)/10“®. 

In Fig. 1 we have plotted the curves up to a redshift of 5. There is another issue to consider, 
however. Our expresions for nc are only valid as long as M' < M 2 (t), due to our approximation 
of loops as point sources. For loops with an initial velocity of vt = 0.25, this condition is satisfied 
only up to a largest redshift of 


50 ’ 


(15) 






(16) 


which is equal to 

^max — (1"^) 

for Vi = 0.25 and Vgq = 0.1 . Beyond that redshift, the values of ncit) and fi„i are suppressed 
beyond the results plotted in Fig. 1 and Eq. since only the tail of the loop ensemble with 
velocities smaller than the mean velocity Vi = 0.25 manage to accrete a substantial amount of 
mass. B is interesting to note that, as mentioned in the introduction, Storrie-Lombardi et al. 
(1996)U find a flattening of Og at z ^ 2 and a possible turnover at z ^ 3. 
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5 Discussion 


We have studied the accretion of hot dark matter onto moving cosmic string loops and made 
use of the results to study early structure formation in the cosmic string plus HDM model. The 
loop accretion mechanism is able to generate nonlinear objects which could serve as the hosts of 
high redshift quasars much earlier than the time cosmic string wakes start becoming nonlinear 
(which for G/i = 10“® and Us = 1/2 occurs at a redshift of about 1). 

The fraction 0„;(z) of the total mass accreted into nonlinear objects by string loops unfor¬ 
tunately depends very sensitively on a and v. On the other hand, this is not surprising since 
the power of the loop accretion mechanism depends on the number and initial sizes of the loops, 
and the scaling relation ~ ua is what should be expected from physical considerations. 

For the values v = \ and a _2 = 1 which are indicated by recent cosmic string evolution 
simulationsEB, we conclude that the loop accretion mechanism produces enough large mass 
protogalaxies to explain the observed abundance of z < 4 quasars (see Figure 1). Note that the 
amplitude of the predicted protogalaxy density curves depends sensitively on the parameters of 
the cosmic string scaling solution which are still poorly determined. Hence, the important result 
is that there are parameters for which the theory predicts a sufficient number of protogalaxies. 
Since not all protogalaxies will actually host quasars, and since the string parameters are still 
uncertain, it would be wrong to demand that the amplitude of the no curve agree with that of 
the observed ng; rather, it should lie above the ng curve. 

It is more difficult to make definite conclusions regarding the abundance of damped Lyman 
alpha absorption systems. In the form of Eq. 15, the condition for the cosmic string loop 
accretion mechanism to be able to explain the data is also satisfied. However, Eq. refers to 
the value of H in baryonic matter. The corresponding constraint on the total matter collapsed 
in structures associated with damped Lyman alpha systems is 


Hdl(^ < 3) > 4 -'10 


where ft is the local fraction of the mass in baryons. From Eq. 5 it follows that the above 
constraint is only marginally satisfied, and this only if the local baryon fraction ft exceeds the 
average value for the whole Universe of about 4 = 0.1. But in the cosmic string model with 
HDM we might expect 4 in nonlinear objects to be enhanced over the average 4 because after 
teq baryons are able to cluster during the time that the HDM is prevented from accreting by the 
free streaming. Thus, cosmic strings may be able to restore agreement with ( |^) in a natural 
way. More calculations are required to resolve this issue. 

Here we have only reported on the mechanism of forming early nonlinear objects through 
accretion onto string loops. Another mechanism is through small-scale structure on the long 
strings leadi^ to the formation of filaments rather than wakes, which has recently been 
investigated^^. It was found that this could be the most effective mechanism, and for the 
maximal possible amount of small-scale structure, flni ~ 1 can be reached already at a redshift 
of 5. It is clearly important to determine the amount of small-scale structure present on strings. 
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Figure 1: Comparison of the number density of host galaxies marks) inferred from 

the observed number density of quasars brigher than Mb = —26 from the Schmidt et al. 
(1991,1995) survey(‘x’ marks) with the number density no of protogalaxies of mass greater 
than lO^^M© predicted in the cosmic string theory with HDM, for the parameters discussed in 
the text, and = 1. The horizontal axis is the redshift. 




